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Abstract. Chemical tagging of stellar debris from disrupted open clusters and as¬ 
sociations underpins the science cases for next-generation multi-object spectroscopic 
surveys. As part of the Galactic Archaeology project TraCD (Tracking Cluster De¬ 
bris), a preliminary attempt at reconstructing the birth clouds of now phase-mixed thin 
disk debris is undertaken using a parametric minimum spanning tree (MST) approach. 
Empirically-motivated chemical abundance pattern uncertainties (for a 10-dimensional 
chemistry-space) are applied to NBODY6-realised stellar associations dissolved into a 
background sea of field stars, all evolving in a Milky Way potential. We demonstrate 
that significant population reconstruction degeneracies appear when the abundance un¬ 
certainties approach ~0.1 dex and the parameterised MST approach is employed; more 
sophisticated methodologies will be required to ameliorate these degeneracies. 


1. Introduction 

The underlying premise of Galactic Archaeology is that surveys provide a fossil record 
of the evolution of the Milky Way. Mining this record entails the search for sub¬ 
clustering in multi-dimensional (spatial, kinematic, chemical) datasets. For systems 
with long dynamical times, relatively few dimensions are needed to identify clustering; 
e.g., energy-angular momentum phase space alone can identify the building blocks of 
the stellar halo (Brook et al. 2003). Unfortunately, the dominant baryonic component 
of the Galaxy - the thin disk - does not fall into this somewhat ‘straightforward’ regime. 

Our thin disk is thought to have been built by many generations of now-disrupted 
stellar associations, the debris from which having been subsequently scattered/migrated 
by a convolution of processes, including systematic spiral arm- and bar-driven ‘churn¬ 
ing’, and random diffusion-like kinematic heating from giant molecular clouds. Before 
an association has fully disrupted, identifying stellar siblings - i.e., the parent birth 
cloud/association - is relatively straightforward. Spatial, kinematic, and/or phase space 
coherency can be maintained for a few ~100 Myrs (depending upon cluster mass, con¬ 
centration, and galactocentric radius/orbit). Unfortunately, on the ~10 Gyr timescale of 
the thin disk, the combined effects of diffusive scattering and radial migration quickly 
wash out this coherency, making sub-clustering analysis in low-order spatial and kine¬ 
matic dimensions a fruitless endeavour (at least for reconstructing the birth locations of 
the sea of Galactic field stars and searching for our own Sun’s siblings). 

To combat this dimensionality ‘problem’. Freeman & Bland-Hawthom (2002) 
proposed the use of 10-20 dimensions of chemistry-space (or ‘C-space’). Dubbed 
‘chemical tagging’, the principle hinges on the presumption that if the gas clouds from 
which the now-dissolved stellar associations formed were chemically homogeneous. 
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even after dissolution and full eonfiguration- and phase-spaee mixing of the debris, the 
parent elouds’ ehemieal ‘fingerprints’ would remain invariant and identifiable. This 
presumpfion has been shown fo hold empirieally, wifh ehemieal homogeneify eon- 
firmed on an elemenf-by-elemenf basis (af fhe ~0.1 dex level, spanning a range of 
nueleosynfhefie proeesses) for >20 assoeiafions (De Silva el al. 2007b, 2009). 

Wifh homogeneify eonfirmed, fhe firsl ‘blind’ ehemieal fagging experimenls were 
eondueled (Milsehang el al. 2013, 2014). High-resolulion speelroseopie dala for field 
slars were analysed, lo allempl a probabilislie approaeh fo idenlifying elusler popula- 
lions lurking in fhe field. Though fhe sludy presented a means lo analyse fhe dalasels 
from nexl-generalion surveys, no means for explieilly delermining elusfer/assoeialion 
reeovery pereenlages was presented. Our goal wilhin TraCD (Traeking Cluster De¬ 
bris: Moyano Loyola el al. 2015) is lo build on Ibis pioneering work and eharaelerise 
paramelrie and non-paramelrie approaehes lo mulfi-dimensional group finding wilhin 
C-spaee, wifh fhe goal being fhe developmenl of fools whieh ean inform upeoming 
surveys sueh as GALAH, WEAVE, and 4MOST. 


2. Method 

As delailed by Moyano Eoyola el al. (2015), our framework is a slalie 3-eomponenl 
(logarilhmie halo, Plummer sphere bulge, Miyamolo-Nagai disk) pofenlial; ~10^ disk 
slars, equally spaeed in ages up lo 10 Gyrs old, are evolved wilh an N-body integra¬ 
tor wilh Irealmenls of bolh random moleeular eloud seallering and syslemalie spiral 
arm ehurning (Sellwood & Binney 2002) applied al eaeh limeslep. This baekground 
sea of slars possess kinemalies eonsislenl wilh Ihose of Ihe Milky Way. We employ 
four NBODY6 realisalions of 250 M© stellar assoeiafions, eaeh evolved in Ihe same 
3-eomponenl polenlial as Ihe baekground slars; Ihese are injeeled into Ihe polenlial al 
various galaeloeenlrie radii ro ranging from 4 to 10 kpe; e.g., Ihe dissolulion lime for a 
elusler injeeled al Ihe solar eirele (8 kpe) is ~0.5 Gyrs. As slars eseape Ihe dissolving 
assoeialion, Iheir Irajeelories are Iraeked wilh Ihe same integrator advaneing Ihe posi- 
lions of Ihe baekground slars, and Ihe same diffusion and ehurning Irealmenls applied. 

We lag our baekground field slars wilh empirieal radial abundanee patterns drawn 
from Euek & Eamberl (2011), using 10 dimensions of C-spaee (Al ,Mg, Si, Ca, Ti, Cr, 
Co, Ni, Y, Nd). Using (of order) Ibis number of elemenls (Vc~10) provides Ihe lever¬ 
age to span Ihe breadlh of nueleosynlhesis sites, while minimising Ihe seareh Ihrough 
parameter spaee. Beeause of Ihe imposed (i) radial abundanee gradienls, and (ii) eon- 
slanl slar formalion history for our disk, Ihe applied diffusion and ehurning means lhal 
we eonsequenlly also impose a temporal evolulion pattern to Ihe abundanee gradienl 
(Gibson el al. 2013). The gradienl’s dependenee on lime to and radial posilion ro is 
used to lag Ihe mean abundanee of Ihe assoeiafions injeeled into Ihe baekground slars al 
any given lime and loealion, wilh Ihe user also imposing an elemenl-by-elemenl seal- 
ler cTc lo eaeh pre-disrupled system. An imporlanl differenee for Ihe latter is lhal Ihe 
ehemieal ‘fingerprinl’ pattern imposed on Ihe pre-disrupled system is nol ‘random’, bul 
instead unique and homogeneous, as per De Silva el al. (2007b, 2009). 

Through use of a minimum spanning free algorilhm (MST - Allison el al. (2009)), 
we allempl to idenlify Ihe debris of disrupted satellites amongsl Ihe baekground stellar 
disk. Erom Ihe 10-dimensional C-spaee, Ihe MST determines a level of similarity Sc 
belween all single stellar eomponenls: 
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where i and j are respective stars and Ac the abundance for element C (Mitschang et al. 
2013). In order to deconvolve the matrix of 6c values into likely parent stellar clusters, 
MST begins building a similarity tree. To do so, stellar components of greatest 6c value 
are joined as nodes, with subsequent iterations joining less similar stellar components 
until all stars are placed in the similarity tree. Having created the tree, a parametric exit 
condition is defined in which the similarity tree is pruned to a value of 6c where the 
number of clusters, k, is present. Fig 1 illustrates the building of a similarity tree for a 
toy distribution, in which it is built (right panel) and subsequently pruned. Due to the 
nature of the similarity matrix, the parameterised pruning method may dissociate tree 
constituents. Such dissociation events thus have no cluster association. 



Figure 1. Toy model illustrating three clustered populations in x - y space (left). 
MST builds a similarity tree, iteratively joining nodes of greatest 6c until all compo¬ 
nents are included in the tree (right). A parametric cluster condition is then applied 
(dashed line) to define the cluster constituents. 


3. Results 

As a proof-of-concept, we evolved four systems with vq (kpc) = [4,6,8,10] for 5 Gyrs. 
To mimic observations, we filter stars at the end of the simulation to only include those 
within 3 kpc of our imposed solar neighbourhood (centred on (x,y, z)=(-8.5,0,0) kpc. 
To determine the success of the MST cluster ‘reconstruction’, four chemical abundance 
uncertainties ctc (dex) were explored: [0.01,0.05,0.10,0.15]. Having tagged each star 
particle (see Fig 2), the MST was deployed, using a parametric exit condition oi k - 

4. From comparisons of the ‘real’ cluster constituents vs the MST associations, the 
successful cluster recovery rates and dissociation percentages were derived (see Tbl 1). 

4. Conclusions 

We employ a parametric form of the MST algorithm to search for now-dissolved stellar 
associations, phase-mixed with a background exponential disk of field stars. Both the 
background sea of stars and the associations themselves were chemically tagged with 
empirically-motivated abundance patterns and a variety of chemical uncertainties and 
association injection radii were explored. Not surprisingly, in the absence of chemi¬ 
cal uncertainties, the association reconstruction accuracy is high; equally unsurprising, 
the accuracy drops dramatically with increasing abundance uncertainties. While not 
meant to be exhaustive, the work presented here is a successful proof-of-concept. Vari¬ 
ous weighting schemes, and non-parametric approaches, urgently need to be explored. 
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Figure 2. Present-day chemistry-radius distributions of four dissolved 250 M© 
stellar associations injected into the disk 5 Gyrs ago at 4 different galactocentric 
radii, with chemical uncertainties of crc = 0.01 dex (left) and ctc = 0.1 dex (right). 


Table 1. MST recovery accuracy and dissociation population percentages for dif¬ 
ferent chemical abundance uncertainties. 


(Tc (dex) Recovery Accuracy (%) Dissociated Population (%) 


0.01 

100 

0 

0.05 

84 

13 

0.10 

50 

17 

0.15 

14 

10 


alongside more sophistieated multi-dimensional group finding algorithms (Sharma & 
Johnston 2012; Mitsehang et al. 2013, 2014). 
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